Corbino-geometry Josephson weak links in thin superconducting films 
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I consider a Corbino-geometry SNS (superconducting-normal-superconducting) Josephson weak 
link in a thin superconducting film, in which current enters at the origin, flows outward, passes 
through an annular Josephson weak link, and leaves radially. In contrast to sandwich-type annular 
Josephson junctions, in which the gauge-invariant phase difference obeys the sine-Gordon equation, 
here the gauge-invariant phase difference obeys an integral equation. I present exact solutions for 
the gauge-invariant phase difference across the weak link when it contains an integral number N of 
Josephson vortices and the current is zero. I then study the dynamics when a current is applied, 
and I derive the effective resistance and the viscous drag coefficient; I compare these results with 
those in sandwich-type junctions. I also calculate the critical current when there is no Josephson 
vortex in the weak link but there is a Pearl vortex nearby. 

PACS numbers: 74.50.+r,74.78.-w,74.25.-q,74.78.Na 
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I. INTRODUCTION 

Thin-film annular Josephson weak links have been 
proposed 1-3 as a test bed for the observation of the in- 
fluence of the Berry phase 4 on the dynamics 5 of a vor- 
tex trapped in the weak link. Recent experiments have 
been carried out by R. H. Hadfield et al. 6 in Corbino- 
geometry thin-film annular Josephson weak links, in 
which the weak links are in the same plane as the elec- 
trodes. The weak links were fabricated using a focused- 
ion-beam technique in a superconductor/normal-metal 
(Nb/Cu) bilayer to mill a 50 nm trench in the super- 
conducting layer to form a weak-link SNS junction. In 
the following, I theoretically examine the properties of 
a thin-film annular Josephson weak link in an idealized 
Corbino geometry, in which current enters at the ori- 
gin, flows outward, passes through an annular Josephson 
weak link, and leaves radially. 

The topological differences between annular weak links 
and straight weak links of finite length produce striking 
differences in behavior. For example, since only inte- 
gral numbers N of flux quanta can be present in annular 
weak links, their critical currents are zero when N ^ 0, 
whereas arbitrary amounts of flux can enter finite-length 
weak links, such that their critical currents are usually 
continous functions of the applied magnetic field. 

I consider here only thin films of thickness d less than 
the London penetration depth A, in which the current 
density j is practically uniform across the thickness, and 
the characteristic length governing the spatial distribu- 
tion of the magnetic field distribution is the Pearl length, 7 



A = 2X 2 /d. 



(1) 



Figure 1 shows the Corbino-geometry SNS Josephson 
weak link considered. Current, supplied to the inner su- 
perconducting (S) film at the origin, flows radially out- 
ward and passes through the annular weak link (N) of in- 
ner and outer radii R- = i? — cZat/2 and R + = R + djf/2, 
where d^ "C R, and continues to flow radially outward 
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FIG. 1: Sketch of the thin-film annular Josephson weak link 
of inner and outer radii R- = R — djv/2 and R+ — R + 
djv/2, considered here. (The weak-link width djv is shown 
greatly exaggerated.) Current enters the film at the origin 
and leaves radially at a distance well beyond R+. Dashes 
show the integration contour used to derive Eq. (5). 



through the outer superconducting (S) film. 

For simplicity, I consider only the case for which A 3> 
R. When a Josephson vortex is trapped in the weak link 
or a Pearl vortex 7 is situated in the vicinity of R, the 
magnetic flux 0o = h/2e carried up through the film is 
spread out over an area of order ir A 2 , so that the corre- 
sponding magnetic flux density is very weak. Although 
we can neglect the magnetic field generated by the vor- 
tex, it is essential to take into account the spatial dis- 
tribution of the current density j or the sheet-current 
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density K = jd. 

In thin-film junctions or weak links 8-10 there is a sec- 
ond important length scale, which characterizes the spa- 
tial variation of the gauge-invariant phase across the 
junction, 11-14 

I = 0o/47r^ o A 2 j c = (j) Q /2-Kp KK c , (2) 

in SI units, where j c (assumed to be independent of 9 in 
Fig. 1) is the maximum Josephson current density that 
can flow radially as a supercurrent through the weak link 
and K c = j c d is the maximum Josephson sheet-current 
density. The case I ^> 2nR corresponds to the small- 
junction limit in straight finite-length junctions, and t <C 
2irR the large-junction limit. 15 

The main goals of this paper are to (a) show that when 
there are N Josephson vortices trapped in the weak link, 
the critical current I c is zero for all values of the ra- 
tio £/2wR, (b) present exact static solutions for the 6 
dependence of the gauge-invariant phase difference for 
arbitrary TV for all values of the ratio £/2irR when the 
applied current is zero, (c) examine the dynamics when 
a current is applied, and (d) show how the critical cur- 
rent density is affected by the presence of a nearby Pearl 
vortex when there is no Josephson vortex trapped in the 
weak link. 

In Sec. II, I derive the basic equation for the gauge- 
invariant phase difference (f)(9) across the weak link 
and note that there are three additive contributions to 
<f>'(9) = d(f>(9)/d9 to be considered. I examine in Sec. Ill 
the contribution due to N flux quanta in the weak link, 
in Sec. IV the contribution due to a Pearl vortex pinned 
nearby, and in Sec. V the contribution due to Joseph- 
son currents. In Sec. VI, I derive the integral equations 
connecting (/>' and sin^. I present exact solutions for 
the gauge-invariant phase difference in a thin-film annu- 
lar Josephson weak link containing a single Josephson 
vortex (N = 1) in Sec. VII and an arbitrary number 
N of Josephson vortices in Sec. VIII, and for all cases 
I work out some consequences for the vortex dynamics 
when a net current / is applied. I calculate in Sec. IX 
the critical current of the annular weak link when there 
is a Pearl vortex nearby, and I briefly summarize all re- 
sults in Sec. X. Appendix A contains general expressions 
for the vector potential and sheet-current density gener- 
ated by N flux quanta in a narrow circular slot of radius 
R, Appendix B contains details of the Josephson-current- 
generated sheet current, and Appendix C presents a brief 
comparison with the properties of sandwich-type annular 
junctions. 

II. GAUGE-INVARIANT PHASE DIFFERENCE 

In the context of the Ginzburg-Landau (GL) 
theory, 16 ' 17 the superconducting order parameter can be 
expressed as * = ^o/e* 7 , where *o is the magnitude of 
the order parameter in a uniform sample, / = |4'|/4'o 
is the reduced order parameter, and 7 is the phase. Let 



us assume that the induced or applied current densities 
are so weak that the suppression of the magnitude of the 
superconducting order parameter is negligible, such that 
/ = 1. For a thin film in which d < A the second GL 
equation (in SI units) can be expressed as 

where K = jd is the sheet-current density, A is the vec- 
tor potential, and B = V x A is the magnetic induction. 

With a sinusoidal current-phase relation, the Joseph- 
son sheet-current density in the radial p direction across 
the weak link is K p (9) = K c sm<f>(9), where K c is the 
maximum Josephson sheet-current density and 4>(9) is 
the gauge-invariant phase difference between the inner 
(p < R-) and outer (p > R + ) superconducting banks, 

2n f R+ 

m=l{R-^)-l{R+,e)-— A p (p,9)dp, (4) 

9o Jr_ 

where R± = R ± dfq/2. A simple relation between 4>(9) 
and the sheet-current densities at p = R- and p = R + 
can be obtained by integrating the vector potential A 
around a loop of width a few coherence lengths larger 
than djv enclosing the weak link with one end of the arc 
at 6' = and the other end at 9' = 9, as shown by 
the dashed contour in Fig. 1. (The weak link may cause 
a proximity-induced suppression of the order parameter 
over a distance of the order of the coherence length £ into 
the superconductor. I assume here that £ <C g?at.) Since 
c?tv A, we can neglect the magnetic flux up through 
the loop. Making use of Eq. (3) for those portions of the 
integration along the inner and outer boundaries of the 
weak link, we obtain 

m = m 

+ ^T^f [R+Ke(R+,e')-R-K e (R_,e')}d9'(5) 
0o Jo 

such that the gauge-invariant phase difference obeys 

^ = HB^[R + K e (R + ,0)-R-Ke{R-,0)]- (6) 
d9 <po 

When a current / enters at the origin and there is nei- 
ther a Josephson vortex trapped in the weak link nor a 
Pearl vortex pinned nearby, the sheet-current K has a 
radial component K p = Ij2-Kp but no azimuthal com- 
ponent (Kg = 0), such that the gauge- invariant phase 
difference (f> is independent of 9. Since diy <C R, the ra- 
dial current of the weak link is I — 2irRK c sin <\> to good 
approximation, and the maximum supercurrent that can 
flow without producing a voltage across the weak link is 
the critical current, I c q = 2ixRK c . 

On the other hand, when flux quanta are trapped in 
the weak link or a Pearl vortex is pinned nearby, az- 
imuthal symmetry is destroyed, the radial component of 
the sheet-current density K p varies as a function of 9, and 
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the azimuthal component of the sheet-current density Kg 
has the property that [R + K e (R +1 9)-R_K e (R_, 9)} ^ 0, 
such that <j) varies with 9 according to Eq. (6). The 
net supercurrent carried through the weak link is I = 
Icos'mcj), where the bar denotes the average over 9, and 
the critical current is given by its maximum value, I c = 
I c0 \sin(j)\ max . 

Although there are nonlinearities associated with the 
properties of Josephson weak links, it is important to 
note that Eq. (6) is a linear equation. Just as the net 
sheet-current density K can be written as a linear sum 
of contributions, so also can d<f>/d9 be written as a lin- 
ear sum of contributions. Since we are interested in the 
behavior when flux quanta are in the annulus, pinned 
vortices are nearby, and radial Josephson currents flow, 
we need to calculate the effects of the linear superpo- 
sition of all three of the corresponding contributions to 
the sheet-current density K and the 9 derivative of the 
gauge-invariant phase difference, 



d<p{9) 
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each of which can be calculated from the corresponding 
discontinuity in pKg(p,9) across the annulus at p = R 
[Eq. (6)]. We next examine each of these three contribu- 
tions in turn. 
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FIG. 2: Streamlines of the sheet current Ki„ generated by a 
Pearl vortex at (a;, y) — (p v , 0) = (7?/2, 0) (black point) inside 
the annulus when 7 = and N = 0, obtained as a contour 
plot of the real part of the complex potential Qi n [Eq. (11)]. 
The bold circle shows the annular weak link, and arrows show 
the current direction. 



III. N FLUX QUANTA IN THE ANNULAR 
WEAK LINK 

Suppose that N flux quanta are trapped in the annular 
weak link in the absence of any nearby Pearl vortices or 
any Josephson currents across the junction. Since we are 
considering only the case that A^> R,we can neglect the 
vector potential term in Eq. (3). For p < 7?_, the phase 
7 = const and K = 0. However, for p > R + , the phase 
winds by multiples of 2tt. When there are TV flux quanta 
in the junction, 7 = —N9, which generates the azimuthal 
sheet-current contribution Kg(p) = Ntfio/irpoAp in the 
region p > R + . From Eq. (6) we obtain 



The above current, phase, and field distributions are 
equivalent to those produced by N Pearl vortices whose 
cores are distributed uniformly around the circle of radius 
7?, such that the total magnetic flux carried up through 
the superconducting film is N<p . See Appendix A for 
details. Choosing an integration contour in the shape of 
a circular sector of radius p > R and central angle 9 in- 
stead of the contour shown in Fig. 1, one can show that 
Eq. (8) is valid for any value of A/7?. 

IV. PEARL VORTEX PINNED NEARBY 

Suppose that a Pearl vortex is pinned at (x, y) — {p v , 0) 
either inside the annular weak link (p v < 7?_) or outside 
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FIG. 3: Streamlines of the sheet current K out generated by a 
Pearl vortex at (a;, y) = (p v , 0) = (27?, 0) (black point) outside 
the annulus when 7 = and N = 0, obtained as a contour 
plot of the real part of the complex potential Q out [Eq. (12)]. 
The bold circle shows the annular weak link, and arrows show 
the current direction. Note that the length scales differ from 
those in Fig. 2. 
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(p v > R+ ) , but no flux quanta are trapped in the annular 
slot nor are there any radial Josephson currents across the 
junction. Since V K = and V x K = 0, with the latter 
equation holding to good approximation because the the 
magnetic field can be neglected when R <C A, the method 
of complex potentials and fields can be used to calculate 
the sheet-current density generated in response to the 
Pearl vortex. In general, the complex potential G(() is 
an analytic function of the complex variable £ = x + iy, 
and the corresponding complex sheet current is 



/C(C) = dG(0/d( = K y (x,y)+iK x (x,y). 



(9) 



Since the radial and azimuthal components of K along 
the unit vectors p = x cos 9 + y sin 9 and 9 = y cos 9 — 
xain9 are K p = (H x x + H y y)/p and Kg = (H y x — 
H x y)/p, we have the relation 

£(C) = £(C)C/P = K e0>> °) + iK p(p> e )> ( 10 ) 



where x = pcosQ, y = psin9, p = \/ x 2 + y 2 , and 9 
ta,n~ 1 (y/x). 

When p v < R_ , the complex potential is 



'0 .foft-Pv 



7T/X A V C - Pi 

InC, p > R+, 



it no A 



(11) 



where ( = x + iy, x = pcos9, y = psm9 7 and pi = R 2 / p v 
corresponds to the radial coordinate of an image vortex. 
Figure 2 shows a contour plot of the real part of Gin ■ The 
contours correspond to streamlines of Ki n . 
When p v > R + , the complex potential is 

Gout(C) = 0, P <R-, 

Figure 3 shows a contour plot of the real part of Gout- 
The contours correspond to streamlines of K out . 

Evaluating K in {Q = dG in (()/d( and /C out (C) = 
dGout(0/d( at p = R- and p = R + and using Eq. (6), 
we find that 



d<p{9) 



where 



d6 



P(p v ,0) = l 



= P(Pv,0), 



\pI-M 



pI + i 



2p v cos 9 



(13) 



(14) 



with p v = p v /R, holds for both p v < R_ and p v > R + . 
When p v = or p v = oo, [d(f>{9) / d9] P = 0. However, as 
p v — > R, [d(j>(9) I ' d9\ p — > 1, the same result as Eq. (8) for 
one flux quantum trapped in the annulus (N — 1). 



V. JOSEPHSON CURRENTS 

Let us next focus on the contribution to the sheet 
current K generated by Josephson currents through the 
junction, ignoring the contributions due to flux quanta in 
the annular weak link or a nearby Pearl vortex. To ob- 
tain the equation determining how 4>(9) varies when the 
radial Josephson current K p {R 1 9) = K c sin 4>{9) varies as 
a function of 9, we start by deriving the Green's function 
for this problem, assuming that the current / entering 
at the origin flows through the weak link with a delta- 
function distribution K 0p (R,9) = (I / R)5(9 - 9'). As in 
Sec. IV, we can use the method of complex potentials. 
The required complex potential is 



So(C) 



±i— In 

2tt 



(.C-C'Y 



(15) 



where ( — x + iy = pe lf> , (' = Re lS , and the upper 
(lower) sign holds for p > R (p < R). The corresponding 
sheet current is 



Kn = 



2np 



p\p 2 -R 2 \± 92pRsm(9 -9') 



p 2 + R 2 -2pRcos(9 -9> 



(16) 



where p = xcos9 + ysin9 and 9 = ycos9 — xsm9. As 
p ->• R, we obtain K 0p (R, 9) = (I/R)6{9 - 9') and 



K oe {R±,e) = ± 



2ttR ± 



cot 



(17) 



The complex potential for a general distribution of ra- 
dial Josephson sheet current K p (9) = K c sin </>(#) can be 
obtained from Eq. (15) by replacing / by K c sin cf)(9')Rd9' 
and integrating over 9': 



G(Q=±i^l sin cj>(9') In 

Z7T 



(C-C) 2 



d6'. (18) 



From this expression we find that the radial and az- 
imuthal components of the sheet current associated with 
the Josephson currents are 



KM 9) = 



K c \p 2 -l 



2np 



sin <j>(6')d6' 



kP 2 + 1 - 2/5cos( 



v J -ir P 2 + 1 - 2p 



w sin(l)(6')sm (0-6')d6' 
cos(9-9') 



(20) 



where p = p/R and the upper (lower) sign holds when 
p > R + (p < R-). The terms involving the azimuthal 
components of the sheet current needed in Eq. (6) are 
given by the principal- value integral, 

R±K e (R± ,9) = ±^£ sin <f>(9') cot (^) d9' , 

(21) 

such that the Josephson-current contribution to Eq. (6) 



is 



#(0)i r r . 

— t-^- = - — -+ sir 



^ino:«':.co: (— r )d<>'. (22) 
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VI. GENERAL EQUATIONS 

Combining the contributions from Eqs. (8), (13), and 
(22), we find that that general equation determining the 
angular dependence of the gauge-invariant phase is 



d(f>(9) 
dff 



= N + P(p v ,6) 
R 



- g -f" sin 4>{6>) cot (?—^-)d#. (23) 

We can invert this integral equation by making use of 
1 r /ff -0\ /ff -6"\ inl 



= 2tt6(6 - 6") - 1. 
The result is 



(24) 



sin (f)(9) = sin0 



where 



sin</> = — / sin <j>(6)d6 = — , (26) 
2?r J_ v I c0 

The term involving TV drops out of Eq. (25) because 

(^V = 0. (27) 



f 



cot 1 



Equation (25) can be converted back to Eq. (23) with the 
help of Eq. (24) and j\ P(p v ,9)d9 = 0. 




FIG. 4: Exact solution (f)(6) [Eq. (29)] vs 6 for a Josephson 
vortex (N = 1) in a thin-film annular weak link for l/R = 
0.1, 0.3, 1, 3, and 10. 



^ 3 



l/R = 0.3 
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FIG. 5: Plot of the derivative of the exact solution (f)(9) = 
d(f>(9)/d9 [Eq. (32)] vs 9 for N = 1 and l/R = 0.3, 1, 3, and 
10. 



VII. N = 1 

A. Exact solution for the static case when I — 

When a Josephson vortex is trapped in the annular 
weak link (N = 1) with no Pearl vortex nearby and the 
current / is zero, the Josephson vortex is stationary, and 
the gauge-invariant phase obeys 



d(f>(9) 
d9 



1 + £ifj inmcot (^V'- (28) 



This equation has an exact solution, corresponding to a 
Josephson vortex centered at 9 = 0, 



(f)(9) = 2 tan" 1 



tan(0/2) 



tan(0i/2) 



+ 7T, 



where 



tan(0i/2) = y/l + (R/l.) 2 - R/l, 



(29) 



(30) 



or, alternatively, 



R _ 1 -tan 2 (0!/2) 
7 ~~ 2tan(0i/2) ' 



(31) 



Note that (f>(— n) = 0, and (j)(n) = 2ir; also tan(0i/2) 1 
when I — > oo, and 9\ — > — > when f — > 0. Figure 4 
shows </>(0) vs 9 for a variety of values of l/R. 
Equation (29) yields 



b>(6) = 



d(j>(6) _ tan(0i/2)[l + tan 2 (0/2)] 
d9 ~ tan 2 (0i/2)+tan 2 (0/2) ' 



(32) 



Note that ^'(0) = cot(0i/2) = + (R/l) 2 + R/l and 
(f>'{±ir) = tan(0i/2) = y/1 + (R/l) 2 - R/l. Figure 5 
shows 4>'(9) vs for a variety of values of l/R. 
From Eq. (29) we also obtain 



sin (f)(9) 



2tan(0 1 /2)tan(0/2) 
~tan 2 (0i/2) +tan 2 (0/2)' 



(33) 
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FIG. 6: Plot of the sine of the exact solution sin 0(0) [Eq. 
(33)] vs for N = 1 and l/R = 0.3, 1, 3, and 10. 



which has its maximum (+1) and minimum (-1) at 

H and = +9\. Defining the angular width 9 core of 
the Joscphson core as the range of values for which 
tt/2 < sin 0(0) < 371-/2, we find for N = 1, 

Ocore = 2#i = 4 tan-Vl + W) 2 - . (34) 

When £/_R — > oo, core = 7r, and when l/R -C 1, CO re ~ 
2£/i?. See Figs. 4, 6, and 9. Note also that sin 0(0) = 
at = — 7r, 0, and 7r. Since sin 0(0) is an odd function 
of 0, sin0 = I/Ico = 0. Figure 6 shows sin 0(0) vs for 
a variety of values of l/R. Figures 4-6 and 9 show that 
the width of the Josephson vortex core increases as l/R 
increases and that as l/R — > oo, the Joscphson vortex 
becomes so spread out that its center can be identified 
only as the place where — ir (or an odd multiple of ir). 

The sheet-current distribution generated by the 
Josephson currents can be calculated from Eqs. (19) and 
(20) using the exact solution for 0(0) given in Eq. (29). 
When l/R = oo, we have simply 0(0) = + n, and the 
sheet current has the simple dipole-like behavior K = 
~K c y for p < i?_ and K = K c (R/p) 2 (y cos 20- £ sin 20) 
for p > R + . For finite values of l/R, K can be evaluated 
analytically as in Appendix B or calculated numerically 
from Eqs. (19) and (20). Streamlines of K, generated 
as contour plots of the real part of Q using Eq. (18), are 
shown in Figs. 7 and 8 for l/R = 30 and l/R = 0.3. 
Recall, however, that the total sheet-current distribution 
for N = 1 is the sum of the Joscphson-current contribu- 
tion shown here and the contribution K = 00 o /7r^ o A/9 
for p > i?, as discussed in Sec. III. 

In sandwich-type annular Josephson junctions [see Ap- 
pendix C], the phase obeys a sine-Gordon equation, 
which involves the sine and the second derivative of the 
phase with respect to the coordinate along the junction. 15 
In thin-film annular junctions discussed here, however, 
the sine of the phase obeys an integral equation, ob- 
tained by partial integration of Eq. (25) with sin = 
and P = 0, 




FIG. 7: Dipole-like streamlines of the Josephson-current con- 
tribution to the circulating sheet-current density K for a 
single flux quantum N = 1 trapped in the weak link when 
l/R = 30 and 1 = 0, obtained as a contour plot of the real 
part of the complex Green's function Q [Eq. (18)] using the 
exact static solution 0(0) [Eq. (29)]. The bold circle shows the 
annular weak link, and the arrows show the current direction. 




sin 0(0') = 



i r 

I^RJ-k 



In 



9' -9 



'(9')d9', (35) 



FIG. 8: Streamlines of the Josephson-current contribution 
to the circulating sheet-current density K for a single flux 
quantum N = 1 trapped in the weak link when l/R = 0.3 
and 1 = 0, obtained as a contour plot of the real part of the 
complex Green's function Q [Eq. (18)] using the exact static 
solution 0(0) [Eq. (29)]. The bold circle shows the annular 
weak link, the Josephson vortex is centered at (x, y) = (7?, 0), 
and the arrows show the current direction. 
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which is analogous to corresponding integral equations 
derived in Rcfs. 8 and 9. The exact solution [Eq. (29)] 
obeys this integral equation, as can be verified by evalu- 
ating Eq. (35) with 



b"(9) = 



tan(0i/2)[l - tan 2 (6»i/2)] tan(0/2) sec 2 (6>/2) 



[tan 2 (0i/2)+tan 2 (0/2)] 2 



(36) 



B. Dynamical behavior when I 



To calculate the critical current [see Sec. II] of a small 
Josephson weak link (which corresponds to the case 
R/£ = 0), one usually can start with a non-current- 
carrying static solution <j> for which sin <j > = 0, add a bias 
phase [3, and then compute sin(0 + /3) = cos <f) sin (3 to 
conclude that the critical current is proportional to the 
average |cos 0| . This procedure remains valid here in the 
limit R/£ = 0, and the result is |cos0| = |cos#| = 0, 
which tells us that the critical current is zero in this 
case. However, this procedure fails for finite values of 
R/l because (j>{6) + (3 is not a solution of Eq. (28). No 
static current-carrying state can be generated from the 
exact solution given in Eq. (29); there is no solution cor- 
responding to a stationary Josephson vortex in the pres- 
ence of a current /. In other words, the critical current 
I c of a thin-film annular Josephson weak link is zero for 
all ratios of R/l. 

As soon as a current I is applied, the gauge-invariant 
phase distribution becomes time-dependent and the weak 
link becomes resistive. The behavior is simplest in the 
limit R/l = 0, for which the voltage measured directly 
across the weak link between p = i?_ and R + is, by 
the Josephson relation, V = (h/2e)d(f>/dt = IR n , where 
R n is the normal-state resistance of the annulus. Since 
the phase cj> slips by 2ir with a frequency v, this oc- 
curs because the straight-line phase distribution, given 
by (f)(9) = 9 + 7r at time t = (similar to the dotted line 
for £/R = 10 in Fig. 4), slides rigidly toward negative 
values of 9 with an angular velocity oj = 2ttv, giving rise 
to a voltage V — hv/2e = 4>qv, where v — (R n /<t>o)I. 

For increasing values of the ratio R/l, the time- 
dependent behavior is more conveniently described in 
terms of Josephson- vortex motion using a quasistatic ap- 
proach. The applied current entering at the origin pro- 
duces a uniform sheet-current density Kj = pKj p — 
pi /2-kR at the annulus. The resulting Lorentz force 19 
Fl = —9Fl = —9Kj p 4>o induces the Josephson vor- 
tex to rotate in a clockwise sense around the annulus. 
When R/l^> 1 (l/R -C 1), the Josephson core becomes 
very compact and the dissipation there becomes quite 
large. As a consequence, for the same current /, the vor- 
tex speed v = 2irRV/4>o and the phase-slip frequency 
v = V/(j)o become smaller than in the opposite limit 
l/R~5> 1, which has the effect of reducing the effective 
resistance of the weak link, R e ff = V/I. This behavior 
is similar to that in sandwich- type annular junctions, as 



discussed in Appendix C. 

To show this quantitatively, we first note that the time 
dependence of all quantities calculated from the exact so- 
lution (j)(9) in Eq. (29) can be obtained to good approxi- 
mation by replacing 9 by 9 + uot. The voltage measured 
directly across the weak link between (p,9) = (R-,9) 
and (R+,9) is V(9,t) = (h/2e)d<j)/dt = fovtf \9 + cot). 
The power delivered to the weak link by the external 
current source is therefore 



I(j)QV 

2tt 



t>'(9 + ut)d6 = IV, 



(37) 



where V, the angular average of the voltage, is equal to 
the time-averaged voltage, < V >= hv/2e = (f>ov. How- 
ever, the power dissipated by the ohmic currents across 



out 



= r 

J —IT 



(9 + ut)/R n ]d9 = (V /RnW 2 , (38) 



where the angular average of <// 2 , obtained from Eq. (32), 
is 

W* = ^ r \4>\0)fd9 = v/l + (R/iy. (39) 

Equating the input power P in to the dissipated power 
Pout, we obtain the effective resistance of the annular 
weak link, 

Reff = V/I = R n /W - Rn/Vl + ( (40) 

and the c orrespondin g phase-slip frequency, v = 

(R n /<p )i/Vi + (R/t) 2 - 

When R/l > 1, such that 0' 2 = R/l to good ap- 
proximation, the Josephson core size (~£) becomes much 
smaller than the circumference of the weak link (2irR), 
and it is then appropriate to think of the Josephson vor- 
tex speed v as being determined by a balance between 
the Lorentz force 19 Fl and a viscous drag force 22 rjv. 
Equating the input power P in = F^v to the dissipated 
power P out = rjv 2 , we obtain the viscous drag coefficient 
(units Ns/m), 



V : 



4:n 2 R n R£' 



(41) 



Note that r\ is inversely proportional to the Josephson 
core size. As discussed in Appendix C, this behavior of r\ 
is similar to that in sandwich-type annular junctions, in 
which 7] is inversely proportional to the Josephson pene- 
tration depth A j. 

The above calculations assume that the maximum 
value of the displacement current density across the weak 
link (e r eo/dN)(dV/dt) is much smaller than the maxi- 
mum Josephson current density j c . This approximation 
is equivalent to the requirement that the vortex speed v 
be much smaller than c, where we find for £/R<^ 1, 



<4V3£d 



N 



V 9e r Ad 



s 1/2 
) C ' 



(42) 
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where e r is the relative dielectric constant in the weak 
link and c is the speed of light in vacuum. Note that c 
is the analog of the Swihart velocity 23 in long sandwich- 
type Josephson junctions. 



VIII. EXACT SOLUTIONS FOR ARBITRARY N 

When N equally spaced Josephson vortices are trapped 
in the annular weak link (N = 1, 2, 3, ...) with no Pearl 
vortex nearby and the current / is zero, the Josephson 
vortex is stationary, and the gauge- invariant phase obeys 



dcf>{9) 
dff 



= N - 



r r . 

f- SI 

2ir£, 



1" sin 0(0') cot {t-JL^dg. (43) 

J — 7T 



An exact solution of this equation, corresponding to one 
Josephson vortex centered at = and the others ar- 
ranged around the annulus with equal angular spacing 
A0 = 2tt/N is 



2 tan -1 



tan(7V0/2) 



tan(N9 N /2) 



(44) 



for —ir/N < 9 < ir/N. For outside this interval in 
the positive (negative) direction, multiples of 2-7T must 
be added to (subtracted from) Eq. (44) to make (p(9) 
continuous with the property that </>(7r) — <j){— n) = 2nN. 
Also 



tan(iV6V2) = y/l + (R/N£) 2 - R/Nt, (45) 
or, alternatively, 



R _ 1 - tan 2 (iV0 w /2) 
7 ~ 2tan(7V0Ar/2) ' 



(46) 



Note that tan{N9 N /2) -> 1 when £ -> oo, and 8 N -> 
t/R -+ when I -> 0. 
Equation (44) yields 



b'{8) 



Ntan(N9 N /2) 



tan 2 (N8 N /2) cos 2 (NO/2) + sin 2 (N0/2) 



■ (47) 



Note that <j>'(0) = Ncot(N6 N /2) = ^/N 2 + (R/£) 2 +R/£ 
and 4>'(±tt/N) = Ntan(N6 N /2) = V 'N 2 + (R/£) 2 ~R/£. 
Equation (44) also yields 

. 2tan(^0Ar/2)sin(7V0/2)cos(Af0/2) 

K ' tan 2 (N8 N /2) cos 2 (N8/2) + s\n 2 {N8/2) 

(48) 

which has a maximum (+1) and a minimum (-1) at 8 = 
— 0jv and 8 = +8 N . Defining the angular width 8 core of 
one of the Josephson cores as the range of 9 values for 
which 7r/2 < sin 0(0) < 37r/2 (modulo 2n) 7 we find for 
arbitrary N, 



8 core = 28 N = (4/JV) tan-^Vl + (.R/Nl) 2 - R/N£\. 

(49) 
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FIG. 9: Angular width of the Josephson core 9 CO re (in units 
of tt) vs l/R for N = 1, 2, 3, and 4 [Eqs. (34) and (49)]. The 
dashed line shows the limiting behavior for Nt/R <C 1. 



When t/R -> oo, 8 core = n/N, and when t/R < 1, 
core w 2£/R. See Fig. 9. 

As in the case N = 1, the critical current of the annular 
weak link is zero for all TV. When a current / is applied, 
the effective resistance can be calculated as in Eqs. (37)- 
(40), except that for arbitrary N we have V = N<j) v, 

Pout = [(M 2 /RnW, and 



b' 2 = N 2 ^/l + {R/N£) 2 , 



(50) 



such that the effective resistance of the annular weak link 
is 

Reff = V/I = R n N 2 /^ = R n /^1 + {R/N£) 2 , (51) 

When R/Nt > 1, such that Ij/ 2 = NR/£ to good 
approximation, the Josephson core size (~ t) becomes 
much smaller than the intervortex spacing (2wR/N). In 
this case the effective resistance of the annular weak link 
containing N Josephson vortices is R e jf = NRi, where 
i?i = R n £/R is the effective resistance when N = 1 in 
this limit [see Eq. (40)]. It is also appropriate in this limit 
to think of the Josephson vortex speed v as being deter- 
mined by a balance between the Lorentz force 19 Fl and 
a viscous drag force 22 r\v. Equating the input power per 
vortex Pi n /N — F^v to the dissipated power per vortex 
Pout/N = ?]v 2 , we obtain exactly the same viscous drag 
coefficient as in Eq. (41). 



IX. CRITICAL CURRENT AFFECTED BY A 
NEARBY PEARL VORTEX 

We next consider the behavior when there is no flux 
quantum in the annular weak link (N = 0) but there is 
a Pearl vortex at (x, y) = (p v ,0) either inside the an- 
nulus (p v < R-) or outside (p v > R+). For simplicity 
let us consider only the case for which £/R is so large 
that we can ignore the effect of the Josephson currents 



9 



on d<p{9)/d9. The equation determining the angular de- 
pendence of the gauge-invariant phase is then simply 



d<p{6) 

de 



P(Pv,6), 



(52) 



where p v — p v /R and P{p v ,0) is given in Eq. (14). Inte- 
gration of Eq. (52) yields the gauge-invariant phase dif- 
ference, 



b v (p v ,9) = 9 - 2tan~ 1 [7 t ,tan(fl/2)], 



where 



Pv 



P„ - 1 



(53) 



(54) 



and the constant of integration is chosen such that 
<fiv(Pvi@) — at 9 = 0, the point on the annulus that 
is closest to the Pearl vortex. 

To calculate the critical current [see Sec. II], we note 
that the net supercurrent carried through the weak link 

= (f> v (p v ,6) + (3, where 



is I = Icosin 4>. Noting that 
/3 is a constant bias phase, also is a solution of Eq. (52), 
we obtain supercurrent-carrying solutions for which I = 
7 c0 cos 4> v sin (3. The critical current is then given by the 
simple result, 



h/l 



cO 



COS( 



4V(1 + Ivf 



or 



Ic/Ico = i-pI, 0<p v <l, 

= 1 - Pv 2 , Pv > 1- 



(55) 



(56) 
(57) 



Note that I c = when p v = 1, which corresponds to the 
case that the Pearl vortex has moved into the annular 
junction. This is equivalent to the state N = 1 discussed 
in Sec. VII. 

Equations (55)-(57) are valid only in the limit R/£ = 0. 
To calculate I c for finite values of Rji would require solv- 
ing Eq. (23) for N — at all p v . While this equation can 
be solved perturbatively for small R/£, the corrections 
to Eqs. (55)-(57) are second order in R/£, such that this 
procedure yields only very small increases in the values 
of Ic/Ico for < p v < 1 and p v > 1. How I c is affected 
for small values of IjR (large R/tj remains unknown. 



X. SUMMARY 

In this paper I have reported a detailed study of the 
properties of a Corbino-geometry annular weak link of 
radius R in a superconducting thin film for which the 
Pearl length 7 A = 2\ 2 /d is much larger than R. I have 
considered separately the contributions due to an inte- 
gral number N of flux quanta trapped in the weak link, 
a Pearl vortex pinned nearby, and the Josephson current 
distribution across the weak link. I derived two equiv- 
alent integral equations describing the gauge-invariant 



phase distribution (f>(9) around the annulus, and I de- 
scribed how these integral equations can be transformed 
into each other. I considered the case of N = 1 with no 
nearby Pearl vortex, first presenting an exact solution for 
4>{9) in the static case when 1 = 0, and then discussing 
the dynamic case for 7 > 0, when the Josephson vortex 
rotates around the annulus at constant angular velocity. 
I then briefly discussed the case of an arbitrary num- 
ber N of equally spaced flux quanta trapped in the weak 
link, again presenting an exact solution for the static case 
when 7 = and discussing the dynamic case when I > 0. 
Finally, I calculated the critical current I c of the weak 
link as a function of the position of a nearby Pearl vor- 
tex and showed that I c = when the Pearl vortex falls 
into the weak link. 

I mentioned in the introduction that thin-film annu- 
lar weak links containing trapped vortices have been 
proposed 1-3 as a place to test for the influence of the 
Berry phase on the vortex dynamics. However, in this 
paper I have assumed that the vortex motion is deter- 
mined only by the principle of conservation of energy: 
the vortex speed was obtained by setting the power sup- 
plied to the weak link equal to the power dissipated via 
ohmic currents. I leave it to other authors to discover 
how this treatment may need to be modified to account 
for the influence of the Berry phase. 
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Appendix A: TV flux quanta in a circular slot 

The vector potential describing the magnetic flux den- 
sity B = V x A when there are N flux quanta trapped 
in a narrow annular slot of radius R in an otherwise thin 
film characterized by the Pearl length A = 2X 2 /d can be 
obtained using a procedure similar to that in Ref. 7 with 
the result A = 9Ag(p, z), where 



Ae(p,z) 



2ir 



J (qR)J 1 (q p) 
1 + qA 



e Tqz dq, (Al) 



and the upper (lower) sign holds when z > (z < 0). 
For r = y/p 2 + z 2 > A, B w ±rN(f>o/2irr 2 , where f = 
(pp + zz)/r. The sheet-current density is K = 6Kq{p), 
where 



Ke(p) 



Nfo 



np A 



Jg(qR)Ji(qp) 
1 + qA 



dq 



, (A2) 
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and 9(x) = 1 when x > 1 and Q(x) — when x < 1. 



Appendix B: Sheet-current density 

The Josephson-current-generated sheet-current den- 
sity can be evaluated analytically from the exact solu- 
tion given in Eq. (29) as follows. The complex current 
density JC(() = K y (x,y) + iK x (x,y) = dQ(()d( can be 
obtained by differentiation of Eq. (18). The correspond- 
ing^(C) =K{Q<i/p = Ke{p,e)+iK p { P ,0) [see Eq. (10)], 
where ( = x + iy = pe l ° , x = pcos9, y — psin9, 
p = \J x 2 + y 2 , and 9 = tan _1 (y/a;), is then 



K p (p,6) 



-K c sm6/p 2 , while for p < 1, K e {p,6) = 



/C(C)=±t 



K C R 



2np 



sin 



(Bl) 



where sin 0(0) is given by Eq. (33), (' = Re l6 and the 
upper (lower) sign holds when p > R (p < R). Changing 
variables to u = tan[(#' — 0)/2), r = tan(#/2), and T£ = 
tan(6*i/2) [Eq. (30)], and employing the definition 



at 



R/e 



i + n l + yiTpW' 



(B2) 



reduces the apparent complexity of the resulting integral, 
whose evaluation yields for p = p/R > 1, 



4X c t,Mt 2 + 1) + p(t 2 -1)] 
' p(l+Te) 2 [(a e -p) 2 + (a e +p) 2 T 2 } 

8K c T t T 

- (l + Te) 2 [(ae-p) 2 + (ae+p) 2 T 2 ]' 
and for p = p/R < 1, 

4K c n[(r 2 -l) + a e p(r 2 + l)} 



Kg(p,9) 
K P (p,9) 



Kg(p,9) 
K p (p, 9) 



+ 



(l + Te ) 2 [(l-a e p) 2 + (l + a e p) 2 T 2 ] 

8K c T iT 

(l + r e ) 2 [(l-a e p) 2 + (l + a e p) 2 T 2 } 



(B3) 
(B4) 

,(B5) 
.(B6) 



Since K — xK x + yK y = pK p + 9 Kg, where p = 
x cos 9 + y sin 9 and 9 = y cos 9 — x sin 9, the above results 
also yield K x and K y via K x = K p cos 9 — Kg sin 9 and 
K y = K p sin 9 + Kg cos 9. 

In the limit as p — > 1 (p — > R + or R-), the above 
results reduce to 



Kg(R ± ,9) 

K P ( P ,9) 



2K c Ti(T 2 -Tt) 

(l+n)(T 2 + r 2 )- 
2K c r e T 



(r! + r 2 Y 



(B7) 
(B8) 



However, Kg(R±,9) can be derived more simply from 
Eq. (23) using N = 1, P = 0, and Eqs. (21), (31), 
and (32), while K p (R, 9) can be derived from K p (R, 9) = 
K c sin (f>(9) and Eq. (33). 

When R/t — 0, ag — 0, and ti = 1, the current pattern 
is dipole-like. For p > 1, Kg(p,9) = K c cos9/p 2 and 



-if c cos and K p (p, 9) = —K c sin 0. 
For other values of R/£, the current pattern is dipolc- 
like only at distances p ^> a>e, where, to good approxi- 
mation, 



K e (p,0) 
K P ( P ,9) 



AK C T£ cos 9 



{l+n) 2 p 2 

4K C T£ sin 9 

~ (l+n) 2 p 2 ' 



(B9) 
(B10) 



Appendix C: Comparison with sandwich- type 
annular junctions 

Numerous experimental studies have been carried out 
in annular Josephson junctions, with some of these 
having the so-called Lyngby geometry. 20 These junc- 
tions can be thought of as long ring-shaped Joseph- 
son junctions sandwiched between a pair of supercon- 
ducting washers. In such junctions, there is only one 
length scale, the Josephson penetration depth 21 Aj = 
(</>o/27T/zojcrfe//) 1 ^ 2 , characterizing the spatial variation 
of both the magnetic field and the non-linear core of 
a Josephson vortex as a function of the coordinate x 
along the length of the junction. Here j c is the maxi- 
mum Josephson supercurrent, d e ff — di + 2A, where di 
is the insulating barrier thickness and A is the London 
penetration depth, and it is assumed that Aj 3> d e ff. 
Starting with solutions of the sine-Gordon equation to de- 
scribe the gauge-invariant phase distribution associated 
with a periodic Josephson- vortex array of period L, Leb- 
wohl and Stephen 22 discussed the Lorentz-force-induced 
motion of the array and calculated the resulting viscous 
drag coefficient per unit length of vortex r/. 

The gauge- invariant phase distribution for a sandwich- 
type annular weak link of radius R and width W, where 
d e ff <C W <§C 2ttR, containing a single Josephson vortex 
(TV = 1) can be obtained from Ref. 22 by simply replacing 
L by 2nR. The solution analogous to Eq. (29) is 



<j){9) = 2 sin" 



r ( Re 













(Cl) 



where kK(k) = ttR/Xj, K(k) is the complete elliptic in- 
tegral of the first kind 24,25 of modulus k and parameter 
k 2 , and sn(u|m)is the Jacobian elliptic function of pa- 
rameter m. In the limits of small and large R/Xj, Eq. 
(Cl) reduces to 

4>{9) = 9 + n, R/Xj < 1, (C2) 
= 2sin" 1 [tanh( J R6»/A./)] +7T, R/Xj > 1. (C3) 

The angular dependence of <f>(9) in a sandwich-type an- 
nular weak link is displayed in Fig. 10. 

The critical current of a sandwich-type annular junc- 
tion containing a single Josephson vortex is zero, and 
for small currents it is a good approximation to assume 
that the phase distribution of Eq. (Cl) rotates around 
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FIG. 10: Plots of the solution <j>(6) [Eqs. (C1)-(C3)] vs 9 for 
a Josephson vortex (N = 1) in a sandwich- type annular weak 
link for \j/R = 0.1, 0.3, 1, 3, and 10. 



the annulus with a speed v = Roj = 2nRv. The effective 
resistance of the junction, calculated as in Sec. VII, is 

R eff = V/I = R n /^, (C4) 
where R n is the normal-state resistance of the annular 



junction, <f)' 2 is the angular average of \d<j){6) / d9} 2 , 

~ 2 - ^« (C5) 
0, (C6) 
(C7) 



irXjk ' 

= l, R/\j 



4i? 

-p, R/Xj » i, 



and -E(fc) is the complete elliptic integral of the second 
kind. 24 - 25 

In the limit Xj <C R, for which the Josephson core 
size (~Aj) is much smaller than the circumference of the 
annulus (2nR), it is appropriate to think of the effective 
resistance as arising from a balance between the Lorentz 
force per unit length of vortex and a viscous drag force 
per unit length. In this limit, the viscous drag coefficient 
per unit length (units Ns/m 2 ) is 22 



7] = 



08 



TT 2 R n RW\j ' 



(C8) 



where W is the width of the annular junction. Note that 
r] is inversely proportional to the Josephson core size. 
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